Controlling spatiotemporal chaos and spiral turbulence in excitable media: A review 
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Excitable media are a generic class of models used to simulate a wide variety of natural systems 
including cardiac tissue. Propagation of excitation waves in this medium results in the formation 
of characteristic patterns such as rotating spiral waves. Instabilities in these structures may lead 
to spatiotemporal chaos through spiral turbulence, which has been linked to clinically diagnosed 
conditions such as cardiac fibrillation. Usual methods for controlling such phenomena involve very 
large amplitude perturbations and have several drawbacks. There have been several recent attempts 
to develop low-amplitude control procedures for spatiotemporal chaos in excitable media which are 
reviewed in this paper. The control schemes have been broadly classified by us into three types: (i) 
global, (ii) non-global spatially-extended and (iii) local, depending on the way the control signal is 
appUed, and we discuss the merits and drawbacks for each. 



PACS numbers: 87.19.Hh, 05.45.Gg, 05.45.Jn, 87.18.Hf 



I. INTRODUCTION 

Excitable media denotes a class of systems that share 
a set of features which make their dynamical behavior 
qualitatively similar. These features include (i) the exis- 
tence of two characteristic dynamical states, comprising 
a stable resting state and a metastable excited state^ (ii) 
a threshold value associated with one of the dynamical 
variables characterising the system, on exceeding which, 
the system switches from the resting state to the excited 
state, and (iii) a recovery period following an excitation, 
during which the response of the system to a supra- 
threshold stimulus is diminished, if not completely ab- 
sent [1]. Natural systems which exhibit such features in- 
clude, in biology, cells such as neurons, cardiac myocytes 
and pancreatic beta cells, all of which are vital to the 
function of a complex living organism. Other examples 
of dynamical phenomena associated with excitable me- 
dia include cAMP waves observed during aggregation of 
slime mold, calcium waves observed in Xenopus oocytes, 
muscle contractions during childbirth in uterine tissue, 
chemical waves observed in the Belusov-Zhabotinsky re- 
action and concentration patterns in CO-oxidation reac- 
tion on Pt(llO) surface. Excitation in such systems is 
observed as the characteristic action potential^ where a 
variable associated with the system (e.g., membrane po- 
tential, in the case of biological cells) increases very fast 
from its resting value to the peak value corresponding 
to the excited state, followed by a slower process during 
which it gradually returns to the resting state. 

The simplest model system capable of exhibiting all 
these features is the generic Fitzhugh-Nagumo set of cou- 
pled differential equations: 

de/ dt = e{l — e){e — h) — dg/dt = e{ke — g)^ (1) 

which, having only two variables, is obviously incapable 
of exhibiting chaos. However, when several such sets are 
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coupled together diffusively to simulate a spatially ex- 
tended media (e.g., a piece of biological tissue made up 
of a large number of cells), the resulting high-dimensional 
dynamical system can display chaotic behavior. The gen- 
esis of this spatiotemporal chaos lies in the distinct prop- 
erty of interacting waves in excitable media, which mu- 
tually annihilate on colliding. This is a result of the fact 
that an excitation wavefront is followed by a region whose 
cells are all in the recovery period, and which, therefore, 
cannot be stimulated by another excitation wavefront, as 
for example when two waves cross each other [2] . Interac- 
tion between such waves result in the creation of spatial 
patterns, referred to variously as reentrant excitations (in 
ID), vortices or spiral waves (in 2D) and scroll waves (in 
3D), which form when an excitation wavefront is broken 
as the wave propagates across partially recovered tissue 
or encounters an inexcitable obstacle Q • The free ends of 
the wavefront gradually curl around to form spiral waves. 
Once formed, such waves become self-sustained sources of 
high-frequency excitation in the medium, and usually can 
only be terminated through external intervention. The 
existence of nonlinear properties of wave propagation in 
several excitable media can lead to complex non-chaotic 
spatiotemporal rhythms. Thus, spiral waves are associ- 
ated with periodic as well as quasiperiodic patterns of 
temporal activity. 

However, in this paper, we shall not be discussing the 
many schemes proposed to terminate single spiral waves, 
but instead, focus on the control of spatiotemporally 
chaotic patterns seen in excitable media (in 2 or 3 dimen- 
sions), that occur when under certain conditions, spiral 
or scroll waves become unstable and break up. Various 
mechanisms of such breakup have been identified [4], in- 
cluding meandering of the spiral focus. If the meandering 
is sufficiently high, the spiral wave can collide with itself 
and break up spontaneously, resulting in the creation of 
multiple smaller spirals (Fig. [Tj). The process continues 
until the spatial extent of the system is spanned by sev- 
eral coexisting spiral waves that activate different regions 
without any degree of coherence. This state of spiral 
turbulence marks the onset of spatiotemporal chaos, as 
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FIG. 1: Onset of spatiotemporal chaos in the 2-dimensional 
Panfilov model of hnear dimension L — 256. The initial con- 
dition is a broken plane wave that is allowed to curl around 
into a spiral wave (left). Meandering of the spiral focus causes 
wavebreaks to occur (centre) that eventually result in spi- 
ral turbulence, with multiple independent sources of high- 
frequency excitation (right). 



indicated by the Lyapunov spectrum and Kaplan- Yorke 
dimension [5i]. 

Controlling spatiotemporal chaos in excitable media 
has certain special features. Unlike other chaotic sys- 
tems, here the response to a control signal is not pro- 
portional to the signal strength because of the existence 
of a threshold. As a result, an excitable system shows 
discontinuous response to control. For instance, regions 
that have not yet recovered from a previous excitation 
or where the control signal is below the threshold, will 
not be affected by the control stimulus at all. Also, the 
focus of control in excitable media is to eliminate all ac- 
tivity rather than to stabilize unstable periodic behav- 
ior. This is because the problem of chaos termination 
has great practical importance in the clinical context, as 
the spatiotemporally chaotic state has been associated 
with the cardiac problem of ventricular fibrillation (VF). 
VF involves incoherent activation of the heart that re- 
sults in the cessation of pumping of blood, and is fatal 
within minutes in the absence of external intervention. 
At present, the only effective treatment is electrical de- 
fibrillation, which involves applying very strong electrical 
shocks across the heart muscles, either externally using 
a defibrillator or internally through implanted devices. 
The principle of operation for such devices is to over- 
whelm the natural cardiac dynamics, so as to drive all 
the different regions of the heart to rest simultaneously, 
at which time the cardiac pacemaker can take over once 
again. Although the exact mechanism by which this is 
achieved is still not completely understood, the danger 
of using such large amplitude control (involving ~ kV 
externally and ~ lOOF internally) is that, not only is 
it excruciatingly painful to the patient, but by causing 
damage to portions of cardiac tissue which subsequently 
result in scars, it can potentially increase the likelihood 
of future arrhythmias, (i.e., abnormalities in the heart's 
natural rhythm). Therefore, devising a low-power con- 
trol method for spatiotemporal chaos in excitable media 
promises a safer treatment for people at risk from poten- 



tially fatal cardiac arrhythmias. 

In this paper, we have discussed most of the recent 
control methods that have been proposed for terminating 
spatiotemporal chaos in excitable media [6] . These meth- 
ods are also often applicable to the related class of sys- 
tems known as oscillatory media, described by complex 
Landau-Ginzburg equation [8], which also exhibit spiral 
waves and spatiotemporal chaos through spiral breakup. 
We have broadly classified all control schemes into three 
types, depending on the nature of application of the con- 
trol signal. If every region of the media is subjected to 
the signal (which, in general, can differ from region to re- 
gion) it is termed as global control] on the other hand, if 
the control signal is applied only at a small, localised re- 
gion from which its effects spread throughout the media, 
this is called local control. Between these two extremes lie 
control schemes where perturbations are applied simul- 
taneously to a number of spatially distant regions. We 
have termed these methods as non- global, spatially ex- 
tended control. While global control may be the easiest 
to understand, involving as it does the principle of syn- 
chronizing the activity of all regions, it is also the most 
difficult to implement in any practical situation. On the 
other hand, local control will be the easiest to implement 
(requiring a single control point) but hardest to achieve. 

In the next section we describe a few of the more com- 
monly used models for studying control of spatiotempo- 
ral chaos in excitable media. Section 3 discusses proposed 
methods of global control, while Section 4 discusses other 
spatially extended schemes. The next section deals with 
local control methods, and we conclude with a brief sec- 
tion containing general discussions about chaos control 
and its implications. 

II. MODELS OF SPATIOTEMPORAL CHAOS 
IN EXCITABLE MEDIA 

The generic Fitzhugh-Nagumo model for excitable me- 
dia (Eq. [1]) exhibits a structure that is common to most 
models used in the papers discussed here. Typically, the 
dynamics is described by a fast variable, e(x, t), and a 
slow variable, ^(x, t), the ratio of timescales being given 
by e. The resulting phase space behavior is shown in 
Fig. [2] (left). For biological cells, the fast variable is often 
associated with the transmembrane potential, while the 
slow (recovery) variable represents an effective membrane 
conductance that replaces the complexity of several dif- 
ferent types of ion channels. For the spatially extended 
system, the fast variable of neighboring cells are coupled 
diffusively. There are several models belonging to this 
general class of excitable media which display breakup of 
spiral waves (in 2D) and scroll waves (in 3D), including 
the one proposed by Panfilov [TO] 

de/dt = V^e- f{e)-g, dg/dt = e{e,g){ke- g). (2) 

Here, /(e) is the function specifying the initiation of the 
action potential and is piecewise linear: /(e) = Cie, 
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FIG. 2: (top) Dynamics in the phase- space of the Fitzhugh- 
Nagumo model, with the resulting time evolution of the action 
potential shown in the inset. The resting state corresponds 
to e = 0,^ = 0. (Bottom) The result of applying a positive 
( "+") or negative ( ") additive perturbation of the same du- 
ration to the e variable: "+" control decreases the threshold 
and makes excitation more likely, while " control decreases 
the duration of the action potential and allows the system 
to recover faster. For the duration of the control signal, the 
e-nullcline shifts upward (downward) for positive (negative) 
perturbation as indicated by the dashed (dash-dotted) curve. 



([3]) only in having dg / dt = f{e) — the functional form 
of /(e) and parameter values being as in Ref. [M]. The 
Aliev-Panfilov model [lH is a modified form of the Pan- 
filov model, that takes into account nonlinear effects such 
as the dependence of the action potential duration on the 
distance of the wavefront to the preceding waveback. It 
has been used for control in Refs. 

alia- 

All the preceding models tend to disregard several com- 
plex features of actual biological cells, e.g., the different 
types of ion channels that allow passage of electrically 
charged ions across the cellular membrane. There ex- 
ists a class of models inspired by the Hodgkin-Huxley 
formulation describing action potential generation in the 
squid giant axon, which explicitly takes such details into 
account. While the simple models described above do re- 
produce generic features of several excitable media seen 
in nature, the more realistic models describe many prop- 
erties of specific systems, e.g., cardiac tissue. The general 
form of such models are described by a partial differential 
equation for the transmembrane potential 



dV_ 
~dt 
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DV^V, 



(4) 



where C is the membrane capacitance density and D is 
the diffusion constant, which, if the medium is isotropic, 
is a scalar. lion is the instantaneous total ionic-current 
density, and different realistic models essentially differ in 
its formulation. For example, in the Luo-Rudy I model 
of guinea pig ventricular cells, lion is assumed to be 
composed of six different ionic current densities, which 
are themselves determined by several time-dependent 
ion-channel gating variables whose time-evolution is gov- 
erned by ordinary differential equations of the form: 



'dt 



(5) 



Here, ^oo = (^^/{a^^P^) is the steady state value of ^ and 
'^^ = a^Xf3^ time constant. The voltage-dependent 

rate constants, and are complicated functions of 
V obtained by fitting experimental data. 



for e < ei, /(e) = —€26 + a, for ei < e < e2, and 
/(e) = Cs{e — 1)^ for e > e2. The physically appropriate 
parameters given in Ref. [10] are ei = 0.0026, 62 = 0.837, 
Ci = 20, C2 = 3, C3 = 15, a = 0.06 and k = 3. The func- 
tion e(e,^) determines the time scale for the dynamics of 
the recovery variable: e(e,^) = ei for e < e2, e(e,^) = €2 
for e > e2, and e(e,^) = 63 for e < ei and g < gi with 
gi = 1.8, ei 1/75, 62 1.0, and 63 0.3 

Simpler variants that also display spiral wave breakup 
in 2D include (i) the Barkley model [11]: 

de/dt = \/^e^e-^e{l-e){e-^-^), dg/dt = e-g, (3) 

the appropriate parameter values being given in Ref. [l2[ , 
and (ii) the Bar-Eiswirth model [13], which differs from 



III. GLOBAL CONTROL 

The first attempt at controlling chaotic activity in ex- 
citable media dates back almost to the beginning of the 
field of chaos control itself, when proportional perturba- 
tion feedback (PPF) control was used to stabilize cardiac 
arrhythmia in a piece of tissue from rabbit heart [2Q]. 
This method applied small electrical stimuli, at intervals 
calculated using a feedback protocol, to stabilize an un- 
stable periodic rhythm. Unlike in the original proposal 
for controlling chaos |2l|, where the location of the sta- 
ble manifold of the desired unstable periodic orbit (UPO) 
was moved using small perturbations, in the PPF method 
it is the state of the system that is moved onto the stable 
manifold. However, it has been later pointed out that 
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PPF does not necessarily require the existence of UPOs 
(and, by extension, deterministic chaos) and can be used 
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even in systems with stochastic dynamics [22|. Later, 
PPF method was used to control atrial fibrillation in hu- 
man heart [23]. However, the effectiveness of such control 
in suppressing spatiotemporal chaos, when applied only 
at a local region, has been questioned, especially as other 
experimental attempts in feedback control have not been 
able to terminate fibrillation by applying control stimuli 
at a single spatial location [7]. 

More successful, at least in numerical simulations, have 
been schemes where control stimuli is applied through- 
out the system. Such global control schemes either apply 
small perturbations to the dynamical variables (e or g) 
or one of the parameters (usually the excitation thresh- 
old). The general scheme involves introducing an exter- 
nal control signal A into the model equations, e.g., in the 
Panfilov model [Eq. 



de/dt 



(6) 



for a control duration r. If ^4 is a small, positive pertur- 
bation, added to the fast variable, the result is an effec- 
tive reduction of the threshold (Fig. [21 bottom), thereby 
making simultaneous excitation of different regions more 
likely. In general, A can be periodic, consisting of a 
sequence of pulses. Fig. [3] shows the results of apply- 
ing a pulse of fixed amplitude but varying durations. 
While in general, increasing the amplitude, or the du- 
ration, increases the likelihood of suppressing spatiotem- 
poral chaos, it is not a simple, monotonic relationship. 
Depending on the initial state at which the control sig- 
nal is applied, even a high amplitude (or long duration) 
control signal may not be able to uniformly excite all re- 
gions simultaneously. As a result, when the control signal 
is withdrawn, the inhomogeneous activation results in a 
few regions becoming active again and restarting the spa- 
tiotemporal chaotic behavior. 

Most global control schemes are variations or modi- 
fications of the above scheme. Osipov and Collins [2^ 
have shown that a low- amplitude signal used to change 
the value of the slow variable at the front and back 
of an excitation wave can result in different wavefront 
and waveback velocities which destabilizes the traveling 
wave, eventually terminating all activity, and, hence, spa- 
tiotemporal chaos. Gray [25| has investigated the termi- 
nation of spiral wave breakup by using both short and 
long-duration pulses applied on the fast variable, in 2D 
and 3D systems. This study concluded that while short 
duration pulses affected only the fast variable, long du- 
ration pulses affected both fast and slow variables and 
that the latter is more efficient, i.e., uses less power, in 
terminating spatiotemporal chaos. The external control 
signal can also be periodic [A = Fsin{u;t)]^ in which case 
the critical amplitude Fc required for terminating activ- 
ity has been found to be a function of the signal frequency 

Other schemes have proposed applying perturbations 
to the parameter controlling the excitation threshold, b. 
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FIG. 3: Global control of the 2-dimensional Panfilov model 
with L — 256 starting from a spatiotemporally chaotic state 
(top left). Pseudo-gray-scale plots of excitability e show the 
result of applying a pulse of amplitude A — 0.833 between 
t = 11 ms and 27.5 ms (top centre) that eventually leads 
to elimination of all activity (top right). Applying the pulse 
between t = 11 ms and 33 ms (bottom left) results in some 
regions becoming active again after the control pulse ends 
(bottom centre) eventually reinitiating spiral waves (bottom 
right). 



Applying a control pulse on this parameter {h = 6/, dur- 
ing duration of control pulse;6 = ^o, otherwise) has been 
shown to cause splitting of an excitation wave into a pair 
of forward and backward moving waves [14]. Splitting of 
a spiral wave causes the two newly created spirals to an- 
nihilate each other on collision. For a spatiotemporally 
chaotic state, a sequence of such pulses may cause termi- 
nation of all excitation, there being an optimal time in- 
terval between pulses that results in fastest control. An- 
other control scheme that also applies perturbation to 
the threshold parameter is the uniform periodic forcing 
method suggested by Alonso et al 0, |2^ for control- 
ling scroll wave turbulence in three-dimensional excitable 
media. Such turbulence results from negative tension be- 
tween scroll wave filaments, i.e., the line joining the phase 
singularities about which the scroll wave rotates. In this 
control method, the threshold is varied in periodic man- 
ner [6 = 6o + hfcos{ijot)] and the result depends on the 
relation between the control frequency uj and the spiral 
rotation frequency. If the former is higher than the lat- 
ter, sufficiently strong forcing is seen to eliminate turbu- 
lence; otherwise, turbulence suppression is not achieved. 
The mechanism underlying termination has been sug- 
gested to be the effective increase of filament tension due 
to rapid forcing, such that, the originally negative ten- 
sion between scroll wave filaments is changed to positive 
tension. This results in expanding scroll wave filaments 
to instead shrink and collapse, eliminating spatiotempo- 
ral chaotic activity. In a variant method, the threshold 
parameter has been perturbed by spatially uncorrelated 
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FIG. 4: Spatiotemporal chaos (top row) and its control (bottom row) in the 2- dimensional Luo-Rudy I model with L = 90 
mm. Pseudo- gray-scale plots of the transmembrane potential V show the evolution of spiral turbulence at times T = 30 ms, 90 
ms, 150 ms and 210 ms. Control is achieved by applying an external current density / = 150/iA/cm^ for r = 2.5 msec over a 
square mesh with each block of linear dimension L/K — 1.35 cm. Within 210 msec of applying control, most of the simulation 
domain has reached a transmembrane potential close to the resting state value; moreover, the entire domain is much below the 
excitation threshold. The corresponding uncontrolled case shows spatiotemporal chaos across the entire domain. 



Gaussian noise, rather than a periodic signal, which also 
results in suppression of scroll wave turbulence [27]. 

As already mentioned, global control, although easy 
to understand, is difficult to achieve in experimental sys- 
tems. A few cases in which such control could be imple- 
mented include the case of eliminating spiral wave pat- 
terns in populations of the Dictyostelium amoebae by 
spraying a fine mist of cAMP onto the agar surface over 
which the amoebae cells grow [28|. Another experimen- 
tal system where global control has been implemented is 
the photosensitive Belusov-Zhabotinsky reaction, where 
a light pulse shining over the entire system is used as 
a control signal [29]. Indeed, conventional defibrillation 
can be thought of as a kind of global control, where a 
large amplitude control signal is used to synchronize the 
phase of activity at all points by either exciting a previ- 
ously unexcited region (advancing the phase) or slowing 
the recovery of an already excited region (delaying the 
phase) [soj. 



IV. NON-GLOBAL SPATIALLY EXTENDED 
CONTROL 

The control methods discussed so far apply control sig- 
nal to all points in the system. As the chaotic activity 
is spatially extended, one may naively expect that any 
control scheme also has to be global. However, we will 
now discuss some schemes that, while being spatially ex- 
tended, do not require the application of control stimuli 
at all points of the system. 



A. Applying control over a mesh 

The control method of Sinha et al [sij involving 
suprathreshold stimulation along a grid of points is based 
on the observation that spatiotemporal chaos in excitable 
media is a long-lived transient that lasts long enough to 
establish a non-equilibrium statistical steady state dis- 
playing spiral turbulence. The lifetime of this transient, 
tl, increases rapidly with linear size of the system, L, 
e.g., increasing from 850 ms to 3200 ms as L increases 
from 100 to 128 in the two-dimensional Panfilov model. 
This accords with the well-known observation that small 
mammals do not get life-threatening VF spontaneously 
whereas large mammals do [32] and has been experimen- 
tally verified by trying to initiate VF in swine ventricular 
tissue while gradually reducing its mass [33]. A related 
observation is that non-conducting boundaries tend to 
absorb spiral excitations, which results in spiral waves 
not lasting for appreciable periods in small systems. 

The essential idea of the control scheme is that a do- 
main can be divided into electrically disconnected regions 
by creating boundaries composed of recovering cells be- 
tween them. These boundaries can be created by trigger- 
ing excitation across a thin strip. For two-dimensional 
media, the simulation domain (of size L x L) is divided 
into smaller blocks by a network of lines with the 
block size {L/ KxL/ K) small enough so that spiral waves 
cannot form. For control in a 3D system, the mesh is used 
only on one of the faces of the simulation box. Control 
is achieved by applying a suprathreshold stimulation via 
the mesh for a duration r. A network of excited and sub- 
sequently recovering cells then divides the simulation do- 
main into square blocks whose length in each direction is 



fixed at a constant value L/K for the duration of control. 
The network effectively simulates non-conducting bound- 
ary conditions (for the block bounded by the mesh) for 
the duration of its recovery period, in so far as it absorbs 
spirals formed inside this block. Note that r need not be 
large at all because the individual blocks into which the 
mesh divides the system (of linear size L / K) are so small 
that they do not sustain long spatiotemporally chaotic 
transients. Nor does which is related to the mesh 
density, have to be very large since the transient lifetime, 
tl^ decreases rapidly with decreasing L. The method has 
been applied to multiple excitable models, including the 
Panfilov and Luo-Rudy models (Fig. [4]). 

An alternative method for controlling spiral tur- 
bulence that also uses a grid of control points has been 
demonstrated for the Aliev-Panfilov model. Two layers 
of excitable media are considered, where the first layer 
represents the two-dimensional excitable media exhibit- 
ing spatiotemporal chaos that is to be controlled, and the 
second layer is a grid structure also made up of excitable 
media. The two layers are coupled using the fast vari- 
able but with asymmetric coupling constants, with ex- 
citation pulses travelling \/D times faster in the second 
layer compared to the first. As the second layer consists 
only of grid lines, it is incapable of exhibiting chaotic 
behavior in the uncoupled state. If the coupling from 
the second layer to the first layer is sufficiently stronger 
than the other way round, the stable dynamics of the 
second layer (manifested as a single rotating spiral) over- 
comes the spiral chaos in the first layer, and drives it to 
an ordered state characterized by mutually synchronized 
spiral waves. 



B. Applying control over an array of points 

An alternative method of spatially extended control is 
to apply perturbations at a series of points arranged in 
a regular array. Rappel et al [34] had proposed using 
such an arrangement for applying a time-delayed feed- 
back control scheme. However, this scheme does not con- 
trol spatiotemporal chaos and is outside the scope of this 
review. 

More recently, the authors [35] have used an array of 
control points to terminate spatiotemporal chaos in the 
Panfilov model. Fig. [5] shows the result of applying a 
spatially non- uniform control scheme, which simulates an 
excitation wave traveling over the system, with the same 
wavefront velocity as in the actual medium. The control 
points are placed distance d apart along a regular array. 
At certain times, the control points at one corner of the 
system is stimulated, followed by the successive stimula- 
tion of the neighboring control points, such that a wave 
of stimulation is seen to move radially away from the site 
of original stimulation. This process is repeated after 
suitable intervals. Note that, simulating a traveling wave 
using the array is found to be more effective at control- 
ling spatiotemporal chaos than the simultaneous activa- 
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FIG. 5: Control of the 2-dimensional Panfilov model (L = 
256) using an array of control points with spacing d = 6 and 
strength of control stimulus A = 2.5. Stimulation is started 
at the top left corner (T = ms) and lasts at each control 
point, as the wave reaches that point, for 17.9 ms. By 200 
ms, the spatiotemporal chaos has disappeared. 



tion of all control points. Using a traveling wave allows 
the control signal to engage all high-frequency sources 
of excitation in the spiral turbulence regime, ultimately 
resulting in complete elimination of chaos. If, however, 
the control had only been applied locally the resulting 
wave could only have interacted with neighboring spi- 
ral waves and the effects of such control would not have 
been felt throughout the system. The efficacy of the 
control scheme depends upon the spacing between the 
control points, as well as the number of simulated travel- 
ing waves. Traveling waves have previously been used in 
Ref. [36] to control spatiotemporal chaos, although in the 
global control context with a spatiotemporally periodic 
signal being applied continuously for a certain duration, 
over the entire system. 



V. LOCAL CONTROL OF SPATIOTEMPORAL 
CHAOS 

We now turn to the possibility of controlling spatiotem- 
poral chaos by applying control at only a small localized 
region of the spatially extended system. Virtually all 
the proposed local control methods use overdrive pac- 
ing^ generating a series of waves with frequency higher 
than any of the existing excitations in the spiral tur- 
bulent state. As low- frequency activity is progressively 
invaded by faster excitation, the waves generated by the 
control stimulation gradually sweep the chaotic activity 
to the system boundary where they are absorbed. Al- 
though we cannot speak of a single frequency source in 
the case of chaos, the relevant timescale is that of the 





FIG. 6: (left) Pacing response diagram for 2D Panfilov model (L = 26) showing relative performance of different waveforms. 
The dash-dotted line represents a sine wave and the solid curve represents a wave of biphasic rectangular pulses, such that 
they have the same total energy. Successful control occurs if the effective frequency lies above the broken line representing the 
effective frequency of chaos (/c), as seen for a larger system (L = 500) at times T = 1000 (center) and T = 3800 (right) time 
units, where the control signal is applied only at the center of the simulation domain. The excitation wavefronts are shown in 
white, black marks the recovered regions ready to be excited, while the shaded regions indicate different stages of recovery. 



spiral waves and is related to the recovery period of the 
medium. Control is manifested as a gradually growing 
region in which the waves generated by the control sig- 
nal dominate, until the region expands to encompass the 
entire system. The time required to achieve termination 
depends on the frequency difference between the control 
stimulation and that of the chaotic activity, with control 
being achieved faster when this difference is greater. 

Stamp et al has looked at the possibility of us- 
ing low- amplitude, high-frequency pacing using a series 
of pulses to terminate spiral turbulence. However, using 
a series of pulses (having various waveform shapes) has 
met with only limited success in suppressing spatiotem- 
poral chaos. By contrast, a periodic stimulation protocol 
[38| has successfully controlled chaos in the 2D Panfilov 
model, as well as other models [39^. The key mecha- 
nism underlying such control is the periodic alternation 
between positive and negative stimulation. A more gen- 
eral control scheme proposed in Ref. [40] uses biphasic 
pacing^ i.e., applying a series of positive and negative 
pulses, that shortens the recovery period around the re- 
gion of control stimulation, and thus allows the genera- 
tion of very high-frequency waves than would have been 
possible using positive stimulation alone. A simple argu- 
ment shows why a negative rectangular pulse decreases 
the recovery period for an excitable system. The stim- 
ulation vertically displaces the e-nullcline and therefore, 
the maximum value of g that can be attained is reduced. 
Consequently, the system will recover faster from the re- 
covery period (Fig. [2l bottom). 

To understand how negative stimulation affects the re- 
sponse behavior of the spatially extended system, we can 
use pacing response diagrams (Fig. [6l left) indicating the 
relation between the control stimulation frequency / and 
the effective frequency /e// , measured by applying a 
series of pulses at one site and then recording the num- 
ber of pulses that reach another site located at a distance 
without being blocked by a region in the recovery period. 
Depending on the relative value of and the recovery 



period, we observe instances of n : m response, i.e., m 
responses evoked by n stimuli. If, for any range of /, the 
corresponding /g// is significantly higher than the effec- 
tive frequency of spatiotemporal chaos, then termination 
of spiral turbulence is possible. As shown in Ref. [4Q| , 
there are indeed ranges of stimulation frequencies that 
give rise to effective frequencies that dominate chaotic 
activity. As a result, the periodic waves emerging from 
the stimulation region gradually impose control over the 
regions exhibiting chaos (Fig. [6]). Note that, there is a 
tradeoff involved here. If /g// is only slightly higher than 
the chaos frequency, control takes too long. On the other 
hand, if it is too high the waves suffer conduction block at 
inhomogeneities produced by chaotic activity which re- 
duces the effective frequency, and therefore, control fails. 

Recently, another local control scheme has been pro- 
posed [41] that periodically perturbs the model parame- 
ter governing the threshold. In fact, it is the local control 
analogof the global control scheme proposed by Alonso 
et al [i2[ discussed in section Hill As in the other meth- 
ods discussed here, the local stimulation generates high- 
frequency waves that propagate into the medium and 
suppress spiral or scroll waves. Unlike the global control 
scheme, bf » bo, so that the threshold can be negative 
for a part of the time. This means that the regions in 
resting state can become spontaneously excited, which 
allow very high-frequency waves to be generated. 



VI. DISCUSSION 

Most of the methods proposed for controlling spa- 
tiotemporal chaos in excitable media involve applying 
perturbations either globally or over a spatially extended 
system of control points covering a significant proportion 
of the entire system. However, in most practical situa- 
tions this may not be a feasible option, either for issues 
of implementation, or because of the high power for the 
control signal that such methods would need. Moreover, 
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if one is using such methods in the chnical context, e.g., 
terminating fibrihation, a local control scheme has the 
advantage that it can be readily implemented with ex- 
isting hardware of the Implantable Cardioverter Defibril- 
lator (ICD). This is a device implanted into patients at 
high risk from fibrillation that monitors the heart rhythm 
and applies electrical treatment when necessary through 
electrodes placed on the heart wall. A low-energy con- 
trol method involving ICDs should therefore aim towards 
achieving control of spatiotemporal chaos by applying 
small perturbations from a few local sources. 

However, the problem with most local control schemes 
proposed so far is that they use very high-frequency 
waves to overdrive chaos. Such waves are themselves 
unstable and may break during propagation, resulting in 
reinitiation of spiral waves after the original chaotic activ- 
ity has been terminated. The problem is compounded by 
the existence of inhomogeneities in real excitable media. 
Recently, Shajahan et al [42] have found complicated de- 
pendence of spatiotemporal chaos on the presence of non- 
conducting regions and other types of inhomogeneities 
in an excitable system. Such inhomogeneities make the 



proposed local control schemes more vulnerable, as it is 
known that high-frequency pacing interacting with, e.g., 
non-conducting obstacles, results in wave breaks and sub- 
sequent genesis of spatiotemporal chaos [43] . 

The search is still on for a control algorithm for ter- 
minating spatiotemporal chaos in excitable media, that 
can be implemented using low power, or, that need be 
applied in only a small, local region of the system, and 
which will yet be robust, capable of terminating spiral 
turbulence without the control stimulation itself break- 
ing up subsequently. The payoffs for coming up with such 
a method are enormous, as the potential benefits include 
an efficient device for cardiac defibrillation. 

Acknowledgements: We would like to thank colleagues 
and collaborators with whom some of the work described 
above has been carried out, especially, Rahul Pandit, 
Ashwin Pande, Avishek Sen, T. K. Shajahan, David J. 
Christini, Ken M. Stein, Johannes Breuer and Antina 
Ghose. We thank IFCPAR and IMSc Complex Systems 
Project (X Plan) for support. 



[1] J. Keener and J. Sneyd, Mathematical Physiology 
(Springer- Verlag, New York, 1998). 

[2] Note that, unlike waves in ordinary diffusive media which 
dissipate as they propagate further, excitation waves are 
self-regenerating. . 

[3] J. Jalife, J.M. Anumonwo, M. Delmar and J.M. Davi- 
denko, Basic Cardiac Electrophysiology for the Clinician 
(Futura Publishing, Armonk, 1999). 

[4] For a discussion of the multiple scenarios of spiral wave 
breakup, see M. Bar, M. Falcke and M. Or-Guil, in: Lec- 
ture Notes in Physics 532, Springer, New York, 1999, 
p. 326; F.H. Fenton, E.M. Gherry, H.M. Hastings and 
S.J. Evans, Ghaos 12, 852 (2002). 

[5] R. Pandit, A. Pande, S. Sinha and A. Sen, Physica A 
306, 211 (2002). 

[6] An earlier review, discussing methods proposed till 2002, 
can be found in Ref. [7]. 

[7] D.J. Gauthier, G.M. Hall, R.A. Oliver, E.G. Dixon- 
Tulloch, PD. Woff and S. Bahar, Ghaos 12, 952 (2002). 

[8] I.S. Aranson and L. Kramer, Rev. Mod. Phys. 74, 99 

(2002) . 

[9] A.V. Panfilov and P. Hogeweg, Phys. Lett. A 176, 295 
(1993). 

[10] A.V. Panfilov, Ghaos 8, 57 (1998). 

[11] D. Barkley, M. Kness and L.S. Tuckerman, Phys. Rev. A 

42, 2489 (1990). 
[12] S. Alonso, F. Sagues, and A.S. Mikhailov, Science 299, 

1722 (2003). 

[13] M. Bar and M. Eiswirth, Phys. Rev. E 48, 1635 (1993). 
[14] M. Woltering and M. Markus, ScienceAsia 28, 43 (2002). 
[15] R.R. Aliev and A.V. Panfilov, Ghaos, Solitons & Fractals 
7, 293 (1996). 

[16] H. Sakaguchi and T. Fujimoto, Phys. Rev. E 67, 067202 

(2003) . 

[17] H. Sakaguchi and Y. Kido, Phys. Rev. E 70, 052901 



(2005). 

[18] A.V. Panfilov, Phys. Rev. E 59, 6251 (1999). 

[19] G.H. Luo and Y. Rudy, Girc. Res. 68, 1501 (1991). 

[20] A. Garfinkel, M.L. Spano, W.L. Ditto and J.N. Weiss, 

Science 257, 1230 (1992). 
[21] E. Ott, G. Grebogi and J.A. Yorke, Phys. Rev. Lett. 64, 

1196 (1990). 

[22] D.J. Ghristini and J.J. Gollins, Phys. Rev. E 52, 5802 
(1995). 

[23] W.L. Ditto, M.L. Spano, V. In, J. Neff, B. Meadows, 
J.J. Langberg, A. Bolmann, A. and K. McTeague, Int. J. 
Bif. Ghaos 52, 593 (2000). 

[24] G.V. Osipov and J.J. Gollins, Phys. Rev. E 60, 54 (1999). 

[25] R.A. Gray, Ghaos 12, 941 (2002). 

[26] S. Alonso, F. Sagues and A.S. Mikhailov, Ghaos 16, 
023124 (2006). 

[27] S. Alonso, J.M. Sancho and F. Sagues, Phys. Rev. E 70, 
067201 (2004). 

[28] K.J. Lee, R.E. Goldstein and E.G. Gox, Phys. Rev. Lett. 

87, 068101 (2001). 
[29] A. P. Munuzuri, V. Perez- Villar and M. Markus, Phys. 

Rev. Lett. 79, 1941 (1997). 
[30] R.A. Gray and N. Ghattipakorn, Proc. Natl. Acad. Sci. 

USA 102, 4672 (2005). 
[31] S. Sinha, A. Pande and R. Pandit, Phys. Rev. Lett. 86, 

3678 (2001). 

[32] A.T. Winfree, When Time Breaks Down (Princeton 

Univ. Press, Princeton, 1987). 
[33] Y.H. Kim, A. Garfinkel, T. Ikeda, T.J. Wu, G.A. Athill, 

J.N. Weiss, H.S. Karagueuzian and P.S. Ghen, J. Glin. 

Invest. 100, 2486 (1997). 
[34] W.J. Rappel, F. Fenton and A. Karma, Phys. Rev. Lett. 

83, 456 (1999). 
[35] S. Sridhar and S. Sinha, forthcoming. 
[36] W. Peng- Ye, X. Ping and Y. Hua-Wei, Ghin. Phys. 12, 



9 



674 (2003). 

[37] A.T. Stamp, G.V. Osipov and J.J. Collins, Chaos 12, 
931 (2002). 

[38] H. Zhang, B. Hu and G. Hu, Phys. Rev. E 68, 026134 
(2003). 

[39] A related case of this control scheme is that proposed 
in G. Yuan, G. Wang and S. Chen, Europhys. Lett. 72, 
908 (2005), where the high-frequency periodic signal is 
applied from the boundaries.. 



[40] J. Breuer and S. Sinha, preprint n\m. CD/ 0406047 (2004). 
[41] H. Zhang, Z. Cao, N.J. Wu, H.P. Ying and G. Hu, Phys. 

Rev. Lett. 94, 188301 (2005). 
[42] T.K. Shajahan, S. Sinha and R. Pandit, Phys. Rev. E 

75, 011929 (2007). 
[43] A.V. Panfilov and J.P. Keener, J. Theor. Biol. 163, 439 

(1993). 



